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Abstract 


The intrinsic nuclear shapes deviating from a sphere not only man- 
ifest themselves in nuclear collective states but also play impor- 
tant roles in determining nuclear potential energy surfaces (PES’s) 
and fission barriers. In order to describe microscopically and self- 
consistently nuclear shapes and PES’s with as many shape degrees 
of freedom as possible included, we developed multidimensionally- 
constrained covariant density functional theories (MDC-CDFTs). In 
MDC-CDFTs, the axial symmetry and the reflection symmetry are 
both broken and all deformations characterized by 8),, with even 
fu are considered. We have used the MDC-CDFTs to study PES’s 
and fission barriers of actinides, the non-axial octupole Y32 correla- 
tions in N = 150 isotones and shapes of hypernuclei. In this Review 
we will give briefly the formalism of MDC-CDFTs and present the 


applications to normal nuclei. 


1. Introduction 


The “shape” provides an intuitive understanding of the 
spatial distribution of matter density in a quantum many- 
body system, such as molecules [1,2], atoms [3], atomic 
nuclei [4-6] and mesons [7]. In atomic nuclei, the occur- 
rence of spontaneous symmetry breaking results in various 
shapes which are associated with different spatial symme- 
tries. If the shape of an atomic nucleus deviates from a 
sphere, we call it a deformed nucleus. The nuclear de- 
formation manifests itself in many observable phenomena, 
including the small and large amplitude collective motions 
with the rotation and the fission as typical examples [8, 9]. 

There are mainly two kinds of ways to parametrize the 
nuclear shape [10,11]. One of them is the two-center or 
two-center-like parametrization [12-16]. The other way is 
of one-center and to make a multipole expansion, 


oo A 
RO, p) = Ro 1+ 5 5 Bp Yan (0, p) , (1) 
A=1 p=— À 


where §),, is the deformation parameter. The way of mul- 
tipole expansion is usually used in mean field calculations. 
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Several typical nuclear shapes are shown schematically in 
Fig. 1. 

Most of the nuclei have shapes similar as a spheroid 
and such shapes are often described by (629 — an ax- 
ial quadrupole deformation parameter. The Nilsson 
perturbed-spheroid parameter €2 [17,18] was also often 
adopted for numerical convenience. It has been predicted 
for quite a long time that in some atomic nuclei, the 
non-axial quadrupole (triaxial) deformation (22 (or y in 
the Hill-Wheeler coordinates) plays an important role. A 
static triaxial shape manifests itself by the wobbling mo- 
tion or chiral doublet bands; both have been extensively 
studied from experimental and theoretical sides [19-25]. 
The existence of multiple chiral doublet (MyD) bands 
in one nucleus was predicted by Meng et al. [26] and 
later multiple chiral doublet bands were observed [27-31]. 
Recently, collective Hamiltonians have been proposed to 
study chiral and wobbling modes [32,33]. It is desirable 
to construct similar Hamiltonians (the collective potential 
and mass parameters) based on covariant density func- 
tionals. In Ref. [34], Bengtsson et al. proposed that one 
of the fingerprints of the triaxiality could be the low-spin 
signature inversion which has been also a hot topic in low 
energy nuclear structure studies [35-39]. Moreover, the 
termination of rotational bands is connected with the de- 
velopment of triaxial shapes in atomic nuclei [40] and the 
triaxiality may also play a crucial role in superheavy nuclei 
[41]. 

The octupole shapes characterized by (39 were pre- 
dicted to be very pronounced in many nuclei, see Ref. [42] 
for a review. Low-lying negative-parity levels, connected 
with the ground state bands via strong E1 transitions, in 
actinides and some rare-earth nuclei are related to reflec- 
tion asymmetric shapes with non-zero 83o [43-52]. The ev- 
idence for a strong octupole deformation has been shown 
in an experimental study of electric octupole transition 
strengths in ??°Rn and ?*4Ra [53] and in !4Ba [54]. It 
was also revealed that the octupole deformation may ex- 
plain the Ra puzzle [55]: The residual proton-neutron in- 
teractions for Ra isotopes show an anomalous enhance- 
ment around N = 135 [56,57]. Recently, Liu et al. iden- 
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Bru = 9 b20 > 0 Boo < 0 Bao > 0 
b22 #0 b30 É 0 b32 F 0 Boo > 0 


Fig. 1: (Color online) Some typical nuclear shapes. In the first 
row (from left to right), a sphere (8), > 0), a prolate spheroid 
(620 > 0), an oblate spheroid (G20 < 0) and a hexadecapole 
shape (G20 > 0 and {40 > 0) are shown. In the second row 
(from left to right), a triaxial ellipsoid (G20 > 0 and B22 Æ 0), a 
reflection asymmetric octupole shape (630 4 0), a tetrahedron 
[G32 Æ 0 and Byun = 0 for all other (Aj)’s] and a largely elon- 
gated and reflection asymmetric shape (G20 >> 0 and 830 Æ 0) 
are shown. Taken from Ref. [68]. 


tified multiple chiral doublet bands in “°Br and observed 
octupole correlations between them [31]. This raises an 
interesting topic on possible existence of nuclear chirality- 
parity quartet bands. The occurrence of the non-axial 
octupole 632 deformation has also been investigated [58- 
67]. More discussions about the 832 deformation will be 
given in Section 4. 

Nuclear deformations corresponding to higher-order 
terms with à > 3 in Eq. (1) are important to different ex- 
tents. The hexadecapole deformation, (49 or €4, has been 
included in deformed mean field potentials since 1960s, 
see, e.g., Ref. [18]. The effects of the Bso or eg deforma- 
tion on the alignment of angular momentum and dynamic 
moment of inertia JC?) in some superheavy nuclei were 
recently revealed [69, 70]. 

These various deformations are also very relevant to 
large amplitude nuclear collective motions, e.g., the fission. 
Since the nuclear fission is discovered [71], the description 
of the fission process has been a difficult and challenging 
task. The fission dynamics are mostly determined by the 
barriers which prohibit the dissolving of the nucleus [72]. 
Since the barrier penetration is involved in the fission pro- 
cess, both the height and the width of the fission barrier 
are very crucial in calculating fission probability [73, 74]. 
In particular, one of the forefronts of modern nuclear sci- 
ence is to explore the existence of charge and mass limits of 
atomic nuclei, i.e., to study superheavy nuclei [75-80]. In 
order to describe the synthesis mechanism of superheavy 
nuclei via heavy ion fusion reactions, it becomes more de- 
sirable to have reliable descriptions for fission barriers for 
superheavy nuclei [73,74, 81-94]. Until now many popu- 
lar nuclear structure models have been employed to cal- 
culate nuclear fission barriers, including the macroscopic- 
microscopic model [14, 81, 95-100], the extended Thomas- 
Fermi plus Strutinsky integral method [101, 102], the non- 
relativistic mean field models (Hartree-Fock or Hartree- 
Fock-Bogoliubov method with the Skyrme force [103-114] 
and the Gogny force [115,116]), and the covariant density 
functional theory [68, 103, 117-133]. These models show 


similar trend of the barrier heights but the quantitative 
disagreement is rather large for some nuclei. 


The axial quadrupole deformation (9 characterizes the 
elongation of a nucleus and is the most relevant shape for 
the nuclear fission. Besides 629 and 649 which is connected 
with the neck formed during fission, many other shapes are 
indispensable for studying properties of the fission barrier 
and for locating the static fission path. Let us take ac- 
tinides as examples. Due to shell effects, a double-humped 
fission barrier appears in actinide nuclei [10]. There may 
also appear a third minimum in the PES’s of light actinide 
nuclei (99, 134-138]. In 1970s it was found in calculations 
with macroscopic-microscopic models that the height of 
inner fission barrier can be reduced by the triaxial distor- 
tion [139-141] and that of the outer one can be lowered 
by the reflection asymmetric (RA) distortion [142]. The 
importance of the triaxial and octupole distortions were 
later also confirmed both in non-relativistic [143] and rel- 
ativistic [121,144] mean field calculations. 


From the above discussions, it is clear that many defor- 
mations are crucial in the ground state, low-lying states, 
potential energy surfaces and fission properties of atomic 
nuclei. It is thus very desirable to develop self-consistent 
theoretical models in which as many deformations as pos- 
sible are considered. This is especially true for studying 
atomic nuclei in unknown regions, including nuclei close 
to proton and neutron drip lines, superheavy nuclei, et al. 


Nowadays, one of the most successful self-consistent 
approaches is the covariant density functional theory 
(CDFT). The CDFT has been used to describe both 
properties of ground states and excited states of nuclei 
in the whole nuclear chart ranging from light to super- 
heavy regions [145-155], including PES’s and fission barri- 
ers of heavy and superheavy nuclei [68, 121-133, 148, 156]. 
In the last years, we have developed multidimensionally- 
constrained covariant density functional theories (MDC- 
CDFTs) in which the reflection symmetry and the axial 
symmetry are both broken [68, 123,130,157]. In MDC- 
CDFTs, all Byy’s with even H, i.e., B20, Boa, B30, Bsa, 
Bao, +++, are included self-consistently. In CDFTs, ei- 
ther the meson exchange (ME) or point-coupling (PC) in- 
teractions can be used. Furthermore, the couplings can 
be non-linear (NL) or density-dependent (DD). All these 
four types of relativistic density functionals, i.e., NL-ME, 
DD-ME, NL-PC and DD-PC, were implemented in MDC- 
CDFTs. The pairing correlations are treated with ei- 
ther the Bogoliubov approach or the BCS method. The 
MDC-CDFTs have been applied to the study of fission 
barriers and PES’s of actinide and trans-actinide nuclei 
[123, 124, 130-132, 158, 159], the non-axial octupole Y32 ef- 
fects in four N = 150 isotones [125], the third minima and 
triple-humped barriers in light actinides [133] and shapes 
of hypernuclei [160, 161]. 


We present a Review of the MDC-CDFTs and the ap- 
plications of MDC-CDFTs on actinide and trans-actinide 
nuclei. In Section 2, we briefly give the formalism of the 
MDC-CDFTs. The study of PES’s and fission barriers of 
actinide nuclei and the third barriers in light actinides will 
be discussed in Section 3. In Section 4, we present results 
of non-axial octupole Y32 effects in N = 150 isotones. A 
summary and some perspectives are given in Section 5. 
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2. Formalism 


In the covariant density functional theory, a nucleus is de- 
scribed as a composite of nucleons which interact either 
through exchanges of mesons and photons or by point- 
couplings [145-155]. The nonlinear coupling terms [162- 
164] or the density-dependent coupling constants [165, 166] 
were introduced in order to give a proper description for 
nuclear matter (e.g., the saturation properties). Accord- 
ingly, there are four kinds of relativistic density function- 
als. In this Section, we take the nonlinear point-couplings 
(NL-PC) covariant density functional theory as an exam- 
ple and present briefly the formalism of MDC-CDFT. The 
time-reversal symmetry is assumed for the nuclei in ques- 
tion. 
The NL-PC Lagrangian reads, 


L=(i7,0" (2) 


where the linear coupling term, the nonlinear coupling 
term, the derivative coupling term, and the Coulomb term 
are given as follows, 


Myy Liin Lui Lack 


Loou; 


1 1 1 
Li = 325P3 + zev PV + 50TSPrs + sOrv Pry: 
la ee ae 272 
La = 3fsPs + 175.5 + qviev! 5 
1 1 1 5 
Laer = 308 [ps]? + z0 [3s pv]? F 9ors [3 prs] 
1 F 
+ 5orvlaverv)’; 
1 l-r: 
Loou = gE Fw F e= Aopv- (3) 


M is the nucleon mass and ag, ay, ars, ATV, BS, YS, YV, 
ôs, Ov, Org and dry are coupling constants corresponding 
to different channels. ps and prg are the isoscalar and 
isovector densities; py and pry are the time-like compo- 
nents of isoscalar and isovector currents. These densities 
and currents read, 


ps = Yy, Prs = YPY, pv = Yh, Pry = YTY. (4) 
With the Slater determinants used as trial wave func- 


tions, applying the mean field approximation and the no- 
sea approximation, one can derive the Dirac equation, 


hrr) = erpr(r), (5) 


where 


h=a-p+8[M+S(r)]+V(r), (6) 
is the single-particle Dirac Hamiltonian which consists of 
two potentials, the scalar potential and the vector poten- 
tial, 


S = asps+arsprs 7+ bsp + ss 

+ dsAps + rsrs: T, (7) 
V = avypy tarvpry -7 + Wwevyev 

+ ôy Apv + ôryv Aprv ET. (8) 


We can solve the RMF equations either in the coordi- 
nate r space or in a complete basis. These equations and 
various extensions were already solved in r space for spher- 
ical nuclei, including the RMF equations [167, 168], the rel- 
ativistic Hartree-Bogoliubov (RHB) equations [169-174], 


the relativistic Hartree-Fock (RHF) equations [175-180], 
and the random phase approximation (RPA) based on the 
RMF and RHF models [181-183]. For deformed nuclei, 
however, it is very difficult to do so because, in addition 
to conventional complications related to two-dimensional 
or three-dimensional spatial lattice techniques [184-187], 
one also has to deal with problems of the variational col- 
lapse ' and fermion doubling in the relativistic frame- 
work. Nevertheless, the RMF equations have been solved 
on a grid in cylindrical coordinates for axially deformed 
nuclei [144,188,189]. Recently, Tanimura et al. devel- 
oped a novel method to solve the RMF equations in three- 
dimensional lattice [190]. 

It is more convenient to solve the RMF equations for de- 
formed nuclei by expanding the single-particle Dirac wave 
functions in terms of a complete basis. Up to now, sev- 
eral bases have been proposed and used widely when solv- 
ing the RMF equations, including the axially deformed 
harmonic oscillator basis [191,192], the three-dimensional 
Cartesian harmonic oscillator basis [193-195], the trans- 
formed harmonic oscillator basis [196-200], the Poschl- 
Teller-Ginocchio basis [201], and the Woods-Saxon basis 
[202-209]. Furthermore, Geng et al. have developed a 
reflection asymmetric relativistic mean field (RAS-RMF) 
model in a basis generated by a two-center harmonic os- 
cillator potential [210]; this model has been used to study 
even-even '46-156Sm in which the important role of the oc- 
tupole deformation on shape phase transitions was found 
[211]. In the MDC-CDFTs, we adopted an axially de- 
formed harmonic oscillator (ADHO) basis [191,192]. The 
axially deformed harmonic oscillator basis consists of the 
eigenstates of the Schrödinger equation, 


h2 
-= + Vg(z,p)| Palro) = Ea®alro), (9) 
where r = (z, p) with p = y x? + y? and 
1 
Va(z, p) = 5M (wpe + wee); (10) 


is the axially deformed harmonic oscillator potential with 
wp (wz) being the oscillator frequencies perpendicular to 
(along) the symmetry axis. The solution of Eq. (9) is 
obtained as 


imp 


=e M 
V 2T X 


A complex number Ca is introduced for convenience. Xs, 
is a two-component spinor. @n,(z) and ne (p) are the 
harmonic oscillator wave functions, 


Palro) = Cadbn. (2) Rni (P) (7). (1) 


1 1 2 _ 2 
TO R (Z) eH 02) 
1 2n! p [m] 
mo = Yetta) 
oO) = FM Gap mall \b, 


1The variational collapse is also a problem one has to solve in 
the basis expansion method, see, e.g., Section IIIA and Fig. 3 of 
Ref. [130] and references therein. 
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Oscillator lengths b, and b, are related to the frequencies, 

= 1//Mw, and b, = 1/,/Mw,. The corresponding 
eigenenergy Ey = w,(2n, + |mi| +1) +w.(nz + 1/2) and 
the major quantum number Na = 2np + [m] + nz. 

These basis states are eigenstates of Jj, with eigenval- 
ues Ka = m+ ms. We define the time-reversal state 
of Palro) as Palro) = TPa(ro) with the time-reversal 
operator T = io, K where K is the complex conjugation 
operator. Ka = —Ka and Ta = Taq where Ta = 1 is 
the parity. The deformation of the basis pasis is defined 


through the relations w, = wo exp (-V5/4rBoasis) and 
Wp = Wo exp (V5/16T Basis), where wo = 1/3 


the frequency of the corresponding spherical oscillator. 
When solving the RMF equations, the Dirac spinor with 
four components is expanded in terms of the complete 
a f? Pa (ro) 


basis {®,(ra)} as, 
Lagi Palro) ) l 


where œ = {n7, np, Mı, Ms} and f% and g& are the ex- 
pansion coefficients. A truncation has to be introduced in 
summations of Eq. (14) in practical calculations. Follow- 
ing Ref. [212], for the upper (large) component of the Dirac 
spinor, the basis states with [n./f,+(2n,+|mi|)/fo] < Nf 
are included in the summation where f, = max(b,/bo, 1) 
and f, = max(b,/bo,1) with bo = 1//Mwo. For the 
expansion of the lower (small) component of the Dirac 
spinor, the truncation of Ng = Ny + 1 is made. 

When neither the spatial reflection symmetry nor the 
axial symmetry is kept, basis states with different {K,7} 
should be included in the expansion in Eq. (14). Thus 
the Hamiltonian matrix becomes much larger than ax- 
ial or reflection symmetric cases. Nevertheless, the sim- 
plex operator S = ie~'% can still be used to keep the 
Hamiltonian matrix to be block-diagonal. S is Hermitian 
with $2 = 1 for a fermionic system. The basis states 
can be labelled with the good quantum number S (the 
eigenvalue of S), $6, = S®, = (-1)K«—!/26,. The 
whole Hilbert space can be divided into two subspaces: 
One subspace with S = +1 consists of states ®, with 
Ka = +1/2, —3/2, +5/2, —7/2, ... and the other sub- 
space with S = —1 consists of the corresponding time- 
reversal states with Ka = —1/2, +3/2, —5/2, +7/2, 
Since the time-reversal symmetry is assumed, we can in- 
clude in the expansion only the basis states with S = +1 
for which we set Ca = 1 when expanding the large compo- 
nent and Cy = i for the small one [cf. Eq. (11)]. We can 
then construct the Hamiltonian matrix with S = +1 and 
diagonalize it to obtain the single-particle energies and 
wave functions. The Hamiltonian matrix with S = —1 
can be obtained by applying the time-reversal operation 
on that with S = +1. 

In the MDC-CDFTs, we make Fourier series expansions 
for the densities in Eq. (4) and the two potentials V(r) and 
S(r), 


(wzw?) is 


vila) = ( (14) 


F(p, 9, 2 (15) 


=> fule,z Fe lino). 


u=- 


As discussed in Refs. [213,214], point-group symmetries 
have been widely used in mean field theories for atomic 


nuclei. In MDC-CDFTs, we assume that the nuclear po- 
tentials and densities are invariant under the following op- 
erations: the reflection with respect to the y-z plane (+), 
the reflection with respect to the x-z plane ($y) and the 
rotation of 180° with respect to the z axis (8), i.e., 


SrG(a, y, 2) = o(-2,y, z), (16 
Sy (x, y, Z) = P(x, —Y, z), (17 
Sd(x,y, z) = p(—z,—y, 2). (18) 


These three operations and the identity Î form the V4 


group ?. Under the V4 symmetry, 
BY = Yra, Syn = 1) Yaa, Oh = (1) Yyy- 
(19) 
Applying these operations on Eq. (1), one gets, 
Bru = Pan = (—1)" Pan. (20) 


Therefore 6), = 0 when p is an odd number. Similarly, 
in Eq. (15), fa = fă = fa and fn = 0 for odd n and the 
expansion given in Eq. (15) is then reduced as 


f(p,¥; zZ)= fo( P, Z ae Yo fal P; Z J cos(2ne), 
(21) 
where 
2r 
fo(p,) = Vor dyf (p,p, 2), (22) 
: r 2 23 
flp) = f defloe,2)eos2ng), (23) 


are all real functions of p and z. For the details about how 
to calculate the matrix elements and how to deal with the 
densities and their derivatives, please refer to Refs. [130, 
157]. 

The pairing correlations are crucial in open shell nu- 
clei. Karatzikos et al. have shown that fission barriers 
are influenced by the pairing force very much [217]. In 
the MDC-RMF models [130], the BCS approach are used 
to treat the pairing correlations. In the MDC-RHB mod- 
els [157], the Bogoliubov transformation was implemented. 
Both in MDC-RMF and in MDC-RHB models, we can use 
the 6-force of zero-range or the separable pairing force of 
finite-range [218-220]. 

To obtain a PES, i.e., Æ = E({8,}), one can perform 
constraint calculations [6]. We have proposed a modified 
linear constraint method in which the Routhian is calcu- 
lated as, 


1 
E' = Exur +) g Ouan: 
Au 
In the (n + 1)th iteration, the variable C,, is determined 
by, 


(24) 


n+1 
Gy = =a + kapu (ay - Bau) ’ 


where Pàp is the desired deformation, ky, is a constant 


(25) 


and el is the value in the nth iteration. 


?The V4 (or V) group is the Klein’s four-group (Vierergruppe in 
German) or Klein’s group [215,216]. It consists of four elements: 
{I,01,02,03} where I is the identity, o 2 = I and ojo; = op when 


tAI Fk. 
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The RMF equations are solved iteratively. After a de- 
sired accuracy is achieved, we can calculate various phys- 
ical quantities. For example, the total binding energy of 
the nucleus reads 


E = Jei (a -p+ BM) yk 
k 


1 1 oo l 2 1 2 
t gesPs t GOuP yr ZOTSPTs 5 Orv Pay 
oe a ee ee 
+ g8sPs + 7IsPs + 7 WV 


1 1 
+ 3osesAps F ov evApy 


1 1 1 
+ 3orsprsAprs + 3orvervAprv + svc} 
(26) 


ae Evair + Bemas 


where Epair is the pairing energy and Ee.m. is the energy 
corresponding to the center of mass correction. The in- 
trinsic multipole moments are calculated from the density 
by 


Qan = f Prov (rr Y 9, (27) 
where Y),,(Q) is the spherical harmonics. The deforma- 
tion parameter àp is obtained from the corresponding 
multipole moment by 


Ar 


Pru = SNR 


Qan» (28) 
where R = rọ A!/® is the radius of the nucleus, the pa- 
rameter ro = 1.2 fm and N represents proton, neutron or 


nucleon numbers. 


3. PES and fission barriers of actinides 


In the study of nuclear fission, the potential energy sur- 
face (PES) is a crucial ingredient. Physical observables, 
including the fragment mass distribution and the half-life, 
are determined to a large extent by the PES of a fissile nu- 
cleus. Actinide nuclei, characterized by a double- or even 
triple-humped fission barriers, have been studied and used 
as the benchmark in many theoretical works 3. There have 
been very extensive studies on fission barriers and PES’s 
for actinides as well as superheavy nuclei, see Ref. [148] 
for a review and Refs. [68, 121-133] for recent progress. 
In particular, in the framework of covariant density func- 
tional theory, Rutz et al. revealed the importance of the 
octupole shape on the second fission barrier in actinides 
[144] and Abusara et al. revealed the crucial role of triax- 
iality on the first fission barrier in actinides [121]. 

We have obtained the PES’s of actinide nuclei in three- 
dimensional deformation space, i.e., Æ = E'(820, G22, 830), 
with the MDC-RMF model [123, 130]. We found that for 
the second fission barriers (the outer one in many cases) 
in these actinides, the triaxial shape is very crucial. Both 
the second barrier and the first one are lowered consider- 
ably due to the triaxial distortion as compared with results 


3Note that fission barriers can not be measured directly in exper- 
iment. One usually assumes a parabolic shape for the fission barrier 
and obtain the height and width of the barrier by fitting available ex- 
perimental fission cross sections [221]. Therefore such barrier heights 
should be treated as “empirical” ones. 


from axially symmetric calculations. With the triaxial de- 
formation considered, a better agreement with empirical 
values of the outer barrier heights in actinide nuclei was 
found. By examining these PES’s, one can learn a lot of in- 
teresting and useful information concerning the shape and 
stiffness of the ground and isomeric states, the fission bar- 
rier heights and the lowering effect of triaxial and octupole 
deformations on fission barriers. Furthermore, by examin- 
ing two-dimensional PES’s and one-dimensional potential 
energy curves (PECs) of some light actinides, the third 
minima and third barriers in these surfaces (curves) were 
located and analyzed in detail [133]. We will present some 
of these results in this Section. 


3.1. 24°Pu: A typical example 
3.1.1. Three-dimensional PES 


The MDC-RMF calculations were carried out for 24°Pu in 
the following deformation lattice: 


e B20 = 0.25-1.70 with a step size of 0.05; 
e B22 = 0-0.25 with a step size of 0.01; 
e 639 = 0-0.50 with a step size of 0.05. 


Totally 26 (for 822) x 11 (for 639) x 30 (for B20) = 8580 
points were calculated. The three-dimensional PES is 
shown in Fig. 2; in each sub-figure the value of (9 is 
fixed and the energy of ?4°Pu is shown as a function of 
B22 and (39. In these calculations, we used the effective 
interaction PC-PK1 [222, 223]. More details can be found 
in Refs. [123, 130]. 

From the first three sub-figures in Fig. 2 with 62o = 0.25, 
0.30 and 0.35, one finds that although the ground state 
with Boo œ~ 0.3 is a bit soft against the reflection asym- 
metric distortion, its equilibrium shape is both reflection 
symmetric and axially symmetric. When it is driven to 
be more elongated, 74°Pu becomes softer against the non- 
axial as well as octupole distortions. From sub-figures 
with Bo9 = 0.40-0.65, one finds that there appears two 
symmetric minima with non-zero (622. Although in these 
sub-figures all minima correspond to $39 = 0, their soft- 
ness against the reflection asymmetric distortion changes 
as 62o increases. The inner fission barrier 74°Pu locates at 
Boo ~ 0.60. In sub-figures with 629 = 0.70-0.95, the sit- 
uation becomes much simpler, where the nucleus is again 
axially symmetric and reflection symmetric. Nevertheless, 
the trend to become softer against the octupole distortion 
can be seen from these sub-figures, which indicates that 
the reflection asymmetric shape becomes more relevant. 

As 829 becomes larger, there appears a superdeformed 
minimum which corresponds to the fission isomer of 24°Pu; 
beyond the second minimum there exists the second sad- 
dle point in the potential energy surface. From sub- 
figures with B29 = 0.85-1.10, we see that, 74°Pu is re- 
flection symmetric around the superdeformed minimum 
with 820 ~ 0.95, but becomes softer against the (39 dis- 
tortion. Starting from {oo = 1.15, there exist two extrema 
with non-zero (39. Although these two extrema locate at 
the 622 = 0 axis, i.e., the triaxial deformation is not im- 
portant, 74°Pu becomes quite soft along the 322 direction. 
With (29 further increasing (1.25 < 69 < 1.35), each 
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Fig. 2: (Color online) Three-dimensional potential energy surface, i.e., E = E (b20, 822, 830), for 74*°Pu obtained from MDC-RMF 
calculations (F is normalized with respect to that of the ground state). In each sub-figure, (20 is fixed as indicated and the 
energy E is displayed as a function of 322 and 830. The contour interval is 1 MeV. Taken from Ref. [130]. 
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Fig. 3: (Color online) Two-dimensional potential energy sur- 
faces of ° Pu [E = E(b20, 630)] from MDC-RMF calculations 
without (a) and with (b) the triaxiality included (£ is normal- 
ized with respect to that of the ground state). The values of 
y are shown in (b) for some points. The dash-dotted line rep- 
resents the static fission path. Full and open circles indicate 
the ground state and the fission isomer. Full and open trian- 
gles display the first and second saddle points. The contour 
interval is 1 MeV. Adapted from Ref. [123]. 


minimum of reflection asymmetry splits into two minima 
with 622 Æ 0. This is the lowering effect of triaxiality on 
the outer fission barrier in 740Pu. Around this barrier, the 
energy is lowered by about 1 MeV by the triaxial distor- 
tion. When (69 > 1.6 (the last three sub-figures in Fig. 2), 
240Py becomes axially symmetric. 

From these studies for 74°Pu, we can draw the follow- 
ing conclusions [123,130]: (1) The shape of 74°Pu in the 
ground state and in the superdeformed isomeric state are 
axially symmetric and reflection symmetric; (2) By ex- 
amining the three-dimensional PES one can learn useful 
information on the stability of 74°Pu against the triaxial 
and octupole distortions: Compared to the ground state, 
the fission isomeric state is much stiffer against both the 
non-axial and reflection asymmetric distortions; (3) While 
around the inner barrier it is triaxial and reflection sym- 
metric, 74°Pu is both triaxial and reflection asymmetric 
around the second saddle point with 820 = 1.3; (4) The 
triaxial deformation appears only around the fission bar- 
riers and its effects are more pronounced around the first 
saddle point. 


3.1.2. Two-dimensional PES 


The three-dimensional PES can be projected onto a two- 
dimensional deformation space in several ways. In Fig. 3 
we show two-dimensional PES [E = E(620, 830)| from 
MDC-RMF calculations. When the triaxial distortion is 
allowed, the energy of ?4°Pu reaches the lowest value at 
each (820, 830) point during the iteration. At some points 
one gets non-zero (22; this means that the nucleus favors 
a triaxial shape. The triaxial shape mainly appears in 


two regions as seen in Fig. 3(b). The first region starts 
from the first saddle point with B25 ~ 0.60 and contin- 
ues almost vertically up to 639 + 1.0: The values of 622 
are around 0.06-0.12 and correspondingly, y ~ 7—11°; the 
triaxial distortion lowers the binding energy by 2 MeV. 
The second region is around the outer fission barrier: The 
B22 © 0.02-0.03 corresponding to y + 1-2°. At the second 
saddle point, the binding energy is lowered by around 1 
MeV due to the triaxial shape. 
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Fig. 4: (Color online) Two-dimensional potential energy sur- 
face of *4°Pu [E = E(ß20, 622)| around the inner barrier (E 
is normalized with respect to that of the ground state). The 
dash-dotted line represents the static fission path. The full star 
denotes the first saddle point. The contour interval is 0.5 MeV. 
Taken from Ref. [131]. 
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Fig. 5: (Color online) Two-dimensional potential energy sur- 
face of 74°Pu [E = E(620, G22)] around the second fission bar- 
rier (E is normalized with respect to that of the ground state). 
The dash-dotted line represents the static fission path. The 
full star denotes the second saddle point. The contour interval 
is 0.1 MeV. Adapted from Ref. [131]. 


Next we focus on the PES’s [E = E(6820, 622)| around 
the two saddle points which are shown in Figs. 4 and 5. 
Starting from its ground state which is axially symmetric, 
the nucleus 74°Pu evolves to the isomeric state through the 
triaxial valley. The inner barrier is located at 629 ~ 0.65 
and {$22 ~ 0.06. The isomeric state keeps an axially sym- 
metric shape. With (629 further increasing, the nucleus 
enters again a valley with triaxiality, then arrives at fis- 
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Fig. 6: (Color online) Potential energy curve of 74°Pu [E = 
E(620)| obtained from MDC-RMF calculations with various 
symmetries imposed (E is normalized with respect to that of 
the ground state): (1) The solid black curve [triaxial (TA) and 
reflection asymmetric (RA)]; (2) The red dashed curve [axial 
symmetric (AS) but RA]; (3) The violet dot-dashed line [AS & 
reflection symmetric (RS)]. The grey square (circle) denotes the 
empirical inner (outer) barrier height. Adapted from Ref. [123]. 


sion configurations. The second fission barrier locates at 
(B20, B22, 830) = (1.21, 0.02, 0.37). 


3.1.3. One-dimensional PEC 


One can further project the three-dimensional or two- 
dimensional PES onto one shape degree of freedom and 
obtain one-dimensional potential energy curves (PECs). 
In Fig. 6 we show the PEC from (29 + —0.2 to the fission 
configuration calculated with various symmetries imposed 
for 74°Pu: the axial symmetry (AS) or triaxial (TA) sym- 
metry together with reflection symmetry (RS) or reflection 
asymmetry (RA). 

In Fig. 6 one can clearly see the importance of the tri- 
axiality around the inner fission barrier as well as the 
crucial role of the axial octupole Y39 shape around the 
outer fission barrier: The triaxial deformation reduces the 
height of the inner fission barrier considerably (more than 
2 MeV), thus leading to a good agreement with the empir- 
ical value; the reflection asymmetric shape becomes very 
important beyond the second minimum and lowers the 
second fission barrier. 

In addition, we have revealed in Ref. [123] for the first 
time that the second fission barrier is also lowered by 
the triaxiality and around 1 MeV is gained in energy for 
the second saddle point when the triaxial distortion is al- 
lowed. This results in a better reproduction of the empiri- 
cal height for the outer fission barrier. Needless to say, this 
feature can be revealed only when the axial and reflection 
symmetries are broken simultaneously. 


3.2. Double-humped PES’s and fission barriers of 
even-even actinides 


3.2.1. One-dimensional PECs around the barriers 


Constraint calculations for a three-dimensional PES are 
very time-consuming. From the three-dimensional cal- 
culations of 74°Pu, we can obtain useful experiences on 
the roles of various nuclear shapes for actinide nuclei. 
Guided by the features found in the one-, two- and three- 


dimensional PES’s of 74°Pu, potential energy curves and 
fission barriers of even-even actinide nuclei have been sys- 
tematically studied in Ref. [130]. 

An actinide nucleus has a reflection symmetric but tri- 
axial shape around the first fission barrier. Therefore one 
can make a one-dimensional constraint calculation with 
the triaxial deformation considered but not the octupole 
shape. In the upper panel of Fig. 7 we present PECs cov- 
ering the ground state and the inner fission barrier for 
some even-even actinide nuclei. Dotted curves represent 
the axially symmetric results and solid curves display re- 
sults from the triaxial MDC-RMF calculations. The em- 
pirical fission barrier heights [221] are shown for compar- 
ison. It can be clearly seen that the triaxial distortion 
lowers the first fission barrier of these even-even actinides 
by about 1 to 4 MeV. The agreement of the calculated 
fission barrier heights with the empirical ones is good for 
most of the nuclei studied here but ?°9:?3?Th and 738:240U, 
Possible reasons for these discrepancies will be discussed 
later. 

It is more complicated to determine the second fis- 
sion barrier height because more shape degrees of free- 
dom become relevant and there are often two or more 
static fission paths. What we have done was the follow- 
ing [123,130]: (1) The axial symmetry is assumed and 
a two-dimensional constraint calculation in the (G29, 830) 
plane is made; (2) From the two-dimensional calculation 
the lowest static fission path 6¢%*S* (29) is approximately 
identified; (3) Along this fission path, we perform MDC- 
RMF calculations with $29 constrained to a given value 
and with both octupole and triaxial deformations allowed. 

ini. — Blowest (Bo) and a small 822 are taken as the initial 
deformation parameters; (4) In this one-dimensional PEC, 
the second saddle point is located and the second barrier 
height is extracted. 

PECs in the region of the isomeric state and the outer 
fission barrier for even-even actinide nuclei are shown in 
the lower panel of Fig. 7. In addition, experimental values 
of the energies for fission isomers [224] are also shown by 
short horizontal lines. One finds that the triaxiality can 
lower the second fission barrier by about 0.5 to 1 MeV for 
most of the actinide nuclei investigated here. The fission 
barrier heights obtained in MDC-RMF calculations with 
the triaxial deformation included are in good agreements 
with the empirical values for these even-even actinides. 


3.2.2. Heights of the first and second fission barri- 
ers 


The heights for inner fission barriers B} calculated with 
the MDC-RMF models are given in Fig. 8 (the upper 
panel). The corresponding empirical values are taken from 
Ref. [221] and shown for comparison. In this figure one 
finds that the triaxial distortion lowers systematically the 
inner barrier heights by 1-4 MeV in these actinide nuclei. 
Generally speaking, the agreement between the calculated 
barrier heights and the empirical ones is fairly good; but 
this is not the case for ?30-232Th and ?38U. The inner bar- 
rier heights calculated from our model are smaller than 
the empirical values by around 1 MeV for 2°°:?3?Th if the 
triaxial deformation is not considered. When the triax- 
ial deformation is included, the calculated heights of inner 
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Fig. 7: (Color online) Potential energy curves [E = E(620)] obtained from MDC-RMF calculations for even-even actinide nuclei 
(E is normalized with respect to that of the ground state). The upper panel shows PECs around the ground states and the 
first fission barriers and the lower panel shows PECs around the fission isomers and the second fission barriers. Dotted curves 
represent the axially symmetric results and solid curves display results from the triaxial calculations. The empirical values for 
fission barrier heights are taken from Ref. [221] and shown by red dots. The green horizontal lines show experimental energies 
for fission isomers [224]. Taken from Ref. [130]. 
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Fig. 8: (Color online) The heights of inner (B$, the upper 
panel) and outer (Bp, the lower panel) fission barriers for some 
even-even actinide nuclei. Open (full) symbols denotes results 
without (with) the triaxial distortion allowed. The empirical 
barrier heights are taken from Ref. [221] and shown by grey 
squares. Taken from Ref. [123]. 


barriers are smaller than the empirical values by around 
2 MeV for these two nuclei. Comparing the two panels in 
Fig. 8, one can find that the inner barrier is lower than 
the outer one for 72°:?32Th. That is, the inner fission bar- 
rier is not the primary one and this may bring in some 
uncertainties when the parameters of inner fission barrier 
are determined [105]. Note that the Skyrme-Hartree-Fock- 
Bogoliubov calculations gave similar results for ?30:232Th 
[105] and another calculation with the covariant density 
functional theory also gave a very low inner barrier for 
232Th [121]. For ?3SU, the inner fission barrier height cal- 
culated from the axial calculation is already quite compa- 
rable with the empirical height. The triaxial deformation 
further lowers this barrier by 1.5 MeV and leads to a dis- 
agreement with the empirical value. 


In Fig. 8 we also show the heights of outer fission bar- 
riers BP (the lower panel) and make comparison with em- 
pirical heights. For most of the nuclei we focused on, the 
triaxial distortion lowers the outer fission barrier by 0.5 
to 1 MeV which is about 10 to 20% of the barrier height. 
With the triaxiality considered, the calculation agrees bet- 
ter with the empirical values for most of even-even ac- 
tinides and the only exception is 74%Cm. The reason for 
this discrepancy has been attributed to the competition 
between two or more static fission paths beyond the first 
barrier [123]. 


We have also examined the dependence of the above 
conclusions on the effective interactions used in MDC- 
RMF calculations. It was found that the lowering effect of 
the triaxial shape on the outer fission barrier is still there 
when effective interactions other than PC-PK1 were used. 
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3.3. The third minima and barriers in light ac- 
tinide nuclei 


As discussed in Sections 3.1 and 3.2, the potential energy 
surface of actinide nuclei shows a double-humped feature 
along the fission path. In as early as 1970s, by using 
macroscopic-microscopic models, Möller et al. have pre- 
dicted that a third minimum may exist beyond the second 
barrier in the PES’s of some actinide nuclei [225-227]. The 
existence of such third minima has been used to account 
for the thorium anomaly [226-228]. High resolution fission 
cross section measurements were later performed for 7?7U 
and ?30-233Th and these measurements indicate that there 
exists shallow third minima in these nuclei [229-236]. 


In Ref. [134], the PES’s of Ra to Th isotopes were stud- 
ied by adopting a macroscopic-microscopic model in which 
a modified harmonic oscillator potential was used. A third 
minimum was predicted for many nuclei at very large 
quadrupole deformation and the third potential pocket 
could be as deep as 1.5 MeV. Later Ćwiok et al. inves- 
tigated the potential energy surfaces of some even-even 
actinides and found that quite deep minima may exist in 
many of these nuclei [237]. 


In Refs. [238-242], experimental evidences for the third 
minima and hyperdeformed states in U isotopes were re- 
ported. The excitation energies of the third saddle point 
and the third minimum in these nuclei were deduced as 
Bur ~% 6 MeV and Eyr = 3-4 MeV, respectively. This 
corresponds to third pockets with a depth around 2-3 
MeV. These results support predictions made in Ref. [237], 
but are consistent neither with the experiments concern- 
ing ?30-233Th and 237U [229-235] nor with the theoretical 
results given in Refs. [134, 225-227]. 


The macroscopic-microscopic model with more shape 
degrees of freedom included was used to study PES’s of 
actinide nuclei [137,138]. In particular, the authors of 
Refs. [137, 138] included the 619 deformation which lowers 
the third barrier substantially and leads to the disappear- 
ance of the third pockets in many actinide nuclei. Third 
minima, which are quite shallow and the depth are less 
than 400 keV, were predicted only in ?30:232Th [137, 138]. 
Furthermore, Ichikawa et al. carried out a systematic 
study of PES’s of actinide nuclei by using the finite-range 
liquid-drop model [81] and found that the third pockets 
for light Th and U isotopes are also quite shallow and the 
depths are less than 1 MeV [2483]. 


Besides various versions of macroscopic-microscopic 
models, self-consistent mean field approaches, e.g., 
Hartree-Fock or Hartree-Fock-Bogoliubov models with 
Skyrme forces [104,105,244] and the Gogny force [245— 
247] and the relativistic mean field model [144], have also 
been adopted to investigate PES’s of actinide nuclei. Gen- 
erally speaking, no very deep third minima were predicted 
in such studies. 


The MDC-RMF model has been used to examine the 
occurrence and features of the third minima and the third 
barriers in the PES’s of some light actinides in Ref. [133]. 
We will present these results in this Section. 
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Fig. 9: (Color online) Potential energy surfaces of U and Th isotopes [E = E(620, 830)] calculated by using the MDC-RMF 
models with the relativistic functionals PC-PK1 (a) and DD-ME2 (b) (E is normalized with respect to that of the ground 
state). The contour interval is 0.5 MeV. Taken from Ref. [133]. 
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3.3.1. Two-dimensional PES’s beyond the second 
fission barrier 


We investigated even-even Th and U isotopes, namely, 
226,228 ,230,232 Th and 232,234,236,238 7]. In Fig. 9 the PES’s in 
the (B20, B30) plane (820 zx 1.0-3.0 and 630 y 0.3-1.5) are 
displayed for these nuclei. These PES’s are obtained by 
using effective interactions PC-PK1 and DD-ME2 [248], 
respectively. From Fig. 9, we can locate the second barrier 
and if they exist, the third barrier and the third minimum. 

For 23?,234,236,2381] the two functionals PC-PK1 and 
DD-ME2 give different predictions concerning the exis- 
tence of third minima. There are no noticeable third min- 
imum in the PES’s from PC-PK1, as is seen in Fig. 9(a). 
The second saddle points of ?32:234,236,238U locate around 
Boo ~ 1.2-1.3 and 639 + 0.3-0.4. With b20 and 639 further 
increasing, the energy becomes smaller and smaller along 
the lowest fission path. However, the PES’s (629 > 1.5 
and 639 > 0.5) calculated with DD-ME2 are different 
from those with PC-PK1, as shown in Fig. 9(b). With 
the exception of 7°°U, one finds a third minimum in the 
PES’s of other U isotopes, i.e., ?32234.238U, The third 
pocket is quite deep for ?38U, but very shallow for ?32U 
and 734U, Similar results have been obtained by using the 
macroscopic-microscopic model [137, 138, 243, 249] and the 
Skyrme Hartree-Fock-Bogoliubov model [244]. 

Compared to U isotopes, more pronounced third min- 
ima were predicted in Th isotopes with both PC-PK1 and 
DD-ME2 functionals. For ??°Th, the depth of the third 
pocket could be 1-2 MeV. With N increasing, the height 
of the third barrier and the energy of the third minimum 
both decrease and the difference between them becomes 
smaller, i.e., the third well becomes shallower. These ten- 
dencies have also been discussed in Ref. [243]. For 2°°Th, 
there is a shallow third well with a depth less than 1 MeV. 
For ?32Th, the third minimum disappears completely in 
the PES calculated from the PC-PK1 functional while 
there is a third well, though quite shallow, in the PES 
obtained from the effective interaction DD-ME2. 

From Fig. 9 we have concluded that more pronounced 
third minima are predicted with the functional DD-ME2 
than PC-PK1. In addition to the lowest fission path, one 
can also notice other possible paths along which there ap- 
pear several saddle points and shallow minima. These fine 
structures in the PES’s of actinides deserve more detailed 
investigations. 


3.3.2. One-dimensional PECs around the third 
minimum and the third barrier 


In Fig. 9, we can find from the two-dimensional PES’s 
calculated with effective interactions PC-PK1 or DD- 
ME2 that there appear third barriers and third minima 
in 226,228,230,232TH and 252,234,238[], Next we examine for 
these nuclei the one-dimensional potential energy curves 
which are shown in Fig. 10. The excitation energy of the 
third minimum, the heights of the second and third fis- 
sion barriers (relative to the ground state energy), and 
the depth of the third pocket defined as the difference be- 
tween the excitation energy of the third minimum and the 
height of the third fission barrier, i.e., AE = Bm — Er, 
together with empirical parameters taken from Refs. [221, 
231, 235, 239, 241, 242] are tabulated in Table 1. 
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Table 1: Heights of the second and third saddle points, 
By and Bym, and energy of the third minimum Ey (in 
MeV) relative to the ground state for even-even Th and U 
isotopes. AE = Bym — Eru is defined as the depth of the 
third potential well. The empirical values are taken from 
Refs. [221, 231, 235, 239, 241, 242] and denoted by “Emp”. 
Taken from Ref. [133]. 


Nucleus Parameters Bu Fin Bu AE 
220 Th DD-ME2 8.76 7.37 9.31 1.94 
PC-PK1 7.94 5.44 6.73 1.29 
228Th DD-ME2 8.16 6.69 7.82 1.13 
PC-PK1 7.19 4.72 5.50 0.78 
230Th DD-ME2 7.84 5.97 6.60 0.63 
PC-PK1 6.56 4.01 4.45 0.44 
Emp [221 6.80 
Emp [231 5.75 5.55 6.45 0.90 
232Th DD-ME2 7.53 5.42 5.92 0.50 
Emp [221 6.70 
Emp [235 6.50 5.40 5.70 0.30 
2327] DD-ME2 7.25 — = = 
Emp [221 5.40 
Emp [241 4.91 3.20 6.02 2.82 
234] DD-ME2 7.01 — — — 
Emp [221 5.50 
Emp [239 3.1 
2387) DD-ME2 7.70 3.70 4.81 1.11 
Emp [|221 5.50 
Emp [242 56 36 56 20 


For 2*°Th, with effective interactions PC-PK1 and DD- 
ME2, the third minimum was predicted and the excitation 
energy is 5.44 and 7.37 MeV, respectively. With PC-PK1, 
the second barrier is higher than the third one, whereas 
the opposite is true with DD-ME2. The depth of the third 
potential well calculated with PC-PK1 (DD-ME2) is 1.29 
(1.94) MeV. For ??8Th, the third potential well becomes 
shallower and the depth predicted from MDC-RMF cal- 
culations with PC-PK1 (DD-ME2) is 0.78 (1.13) MeV. 
The third minimum is the most pronounced in ?76?28Th 
among the nuclei considered in Ref. [133]. For ?°°Th, 
Emr = 5.55 MeV and Bir = 6.45 MeV have been de- 
duced in Ref. [231]. The MDC-RMF model reproduced 
these values with DD-ME2, but underestimated them by 
around 2 MeV with PC-PK1. The calculation results with 
PC-PK1 and DD-ME2 are quite different for ???Th: There 
does not appear a third minimum in the PES calculated 
with PC-PK1; however, the opposite is true for the case of 
DD-ME2. A shallow third potential well with AF ~ 0.30 
MeV for ?32?Th has been deduced in electron-induced fis- 
sion cross section measurements [235]. The MDC-RMF 
result with DD-ME2, AF = 0.50 MeV, agrees reasonably 
well with this empirical value. In Table 1, one again finds 
that, with the neutron number N increasing, the third 
potential pocket becomes shallower in the Th isotopes. 

For U isotopes, the situation is much simpler than that 
for Th isotopes: Only in ?38SU and with the effective in- 
teraction DD-ME2, the MDC-RMF calculations predict a 
third minimum. As seen in Table 1, the calculated ex- 
citation energy of this minimum is 3.70 MeV and close 
to the empirical value 3.6 MeV [242]. But the calculated 
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Fig. 10: (Color online) Potential energy curves [E = E (820)] 


to that of the ground state). Taken from Ref. [133]. 
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Fig. 11: (Color online) The single-particle level schemes for neutrons (upper panel) and protons (lower panel) of Th near the 
Fermi surface along the static fission path. For B20 < 0.6, only reflection-symmetric deformations are considered; the red (blue) 
curves represent positive (negative) parity states and the Nilsson quantum numbers are shown for some proton levels. When 
B20 > 0.6, the octupole deformation (39 has a non-vanishing value, thus the parity 7 is not a good quantum number and only 
the projection of the single-particle angular momentum on the symmetry axis is shown (the axial symmetry is assumed, see the 
text for more details). The dash-dotted (green) curves denote the Fermi surface and the red (blue) symbols in the lower panel 
are used to guide the eye. The MDC-RMF calculation has been performed with the functional DD-ME2. Taken from Ref. [133]. 
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height of the third barrier is 4.81 MeV which is less than 
the empirical value 5.6 MeV [242]. Thus the MDC-RMF 
model prediction for the depth of the third potential well 
in 795U, AE = 1.11 MeV, is smaller than the correspond- 
ing empirical value 2.0 MeV [242]. 


3.3.3. Single-particle level structure 


The appearance of a hyperdeformed minimum roots in the 
scheme of single-particle levels. This has been discussed 
in Ref. [133]. Figure 11 shows the single-particle orbitals 
around the Fermi surface as functions of the quadrupole 
deformation 829 for ?2°Th. Only the results with the 
DD-ME2 functional are shown. The quadrupole deforma- 
tion of the superdeformed state of ??6Th locates around 
Bao 0.6. For Bọ < 0.6, the reflection symmetry is 
kept, and thus the parity m is conserved. However, when 
B29 > 0.6, the octupole deformation 639 is non-zero and 
the parity is not a good quantum number. Furthermore, 
the triaxial deformation starts to be important around 
the second saddle point and the projection of the single- 
particle angular momentum on the symmetry axis is not 
a good quantum number either. This could leads to a 
quite complicated single-particle structure around the sec- 
ond saddle point. In Ref. [133], we were mainly focused on 
single-particle orbitals beyond the second fission barrier, 
we plotted single-particle levels obtained in axial MDC- 
RMF calculations in Fig. 11. 

The third minimum of 22°Th locates at 62o & 1.5 and 
839 ~ 0.7. From the neutron single-particle levels shown in 
the upper panel of Fig. 11, one can find that around 820 = 
1.5, although the energy gap around the Fermi surface is 
not very large, the density of single-particle levels is quite 
low. As shown in the lower panel of Fig. 11, there is a big 
energy gap around Z = 90 in the region B20 ~ 1.5. Thus it 
has been concluded in Ref. [133] that the existence of the 
hyperdeformed minimum in the PES of ??6Th is mainly 
due to the large proton gap at Z = 90. It has been noted 
that some other single-particle orbitals [dotted with red 
(blue) symbols and labeled with Q, i.e., the projection of 
of the single-particle angular momentum on the symmetry 
axis in the lower panel of Fig. 11] close to the proton Fermi 
surface are also responsible for the formation of this energy 


gap. 


w 
~ 


4. Y32 correlations in N = 150 isotones 


The pure non-axial octupole 832 deformation (i.e., Bàu = 0 
if A Æ 3 or u Æ 2) has a tetrahedral symmetry. The non- 
trivial irreducible representation of the symmetry group 
TP leads to highly degenerate single-particle levels and 
large shell gaps. It has been predicted that the 832 shape 
may appear in the ground states of nuclei with certain 
combination of neutron and proton numbers [60, 63, 67, 
250]. Recently, many theoretical efforts were devoted to 
the study of the tetrahedral shape in atomic nuclei, ei- 
ther from the TP -symmetric single-particle level schemes 
(60, 63, 251, 252] or from different nuclear models includ- 
ing the macroscopic-microscopic model [63, 64, 252, 253], 
the Skyrme Hartree-Fock (SHF) model [61, 62, 64-66, 252- 
254], and the reflection asymmetric shell model (RASM) 
[255,256]. A negative-parity band in 1°°Gd has been sug- 
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Fig. 12: (Color online) Potential energy curves [E = E(832)] 
for (a) 74°Cm, (b) **8Cf, (c) 2°°Fm and (d) ?°?No from MDC- 
RMF calculations (E is normalized with respect to that of the 
ground state). (b) is adapted from Ref. [131]. 
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Fig. 13: (Color online) The single-particle levels near the Fermi 
surface for protons of 745Cf as functions of quadrupole defor- 
mation 6209 (left side) and of 832 with (20 fixed at 0.3 (right 
side). Taken from Ref. [125]. 


gested to be a candidate manifesting tetrahedral symme- 
try [64]. Later, several experiments were made to inves- 
tigate this feature in °°Gd [257, 258]. Although the exis- 
tence of tetrahedral shape has not been confirmed exper- 
imentally, the interests in nuclear tetrahedral symmetry 
do stimulate more and more theoretical studies related to 
point group symmetries in nuclei; the readers are referred 
to Ref. [259] for a recent review on this topic. 


In recent years, the investigation of nuclei with Z ~ 100 
becomes more and more interesting because such stud- 
ies can, on one hand, reveal the structure of these nuclei 
themselves and on the other hand, provide useful infor- 
mation for the structure of superheavy nuclei [70, 260- 
265]. One interesting topic is about the low-lying 27 
states in several even-Z N = 150 isotones [266]. For 
example, in a quasiparticle phonon model with octupole 
correlations included, the low-lying 27 states in the iso- 
tones with N = 150 was explained by the existence of 
the octupole correlations originating from the neutron 
configuration 9/27 [724] @ 5/2*[622] and proton configu- 
rations 9/27 [633] & 5/27 [521] or 7/2+[633] @ 3/27 [521] 
[267]. Since in these two-quasiparticle configurations, the 
projections of the single-particle angular momenta of like- 
quasiparticles on the symmetry axis (K) differ by 2, Y32 
correlations should be more relevant than Y39 correla- 
tions. Indeed, Chen et al., using the reflection asymmetric 
shell model with Y32 correlations considered, have repro- 
duced quite well the low-energy 2~ rotational bands in the 
N = 150 isotones [256]. Furthermore, they also predicted 
that the strong non-axial octupole Y32 effects may still 
exist in superheavy nuclei with Z = 108 [256] and even 
heavier nuclei [268]. Note that in the reflection asymmet- 
ric shell model and many other models, the values of (39 
are input parameters. 

We have carried out a self-consistent and microscopic 
study of Y32 effects in the N = 150 isotones by using the 
MDC-RMF model [125]. It was found that 832 > 0.03 and 
the energy gain due to the inclusion of the 832 deformation 
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E (MeV) 


Fig. 14: (Color online) The single-particle levels near the Fermi 
surface for neutrons of 748Cf as functions of quadrupole defor- 
mation 8209 (left side) and of 832 with (20 fixed at 0.3 (right 
side). Taken from Ref. [125]. 


is more than 300 keV for the ground states of 74°Cf and 
250Fm. In 246Cm and 752No, there appear shallow 832 
minima compared to 748Cf and 75°Fm. The origin of the 
occurrence of the Y32 correlations was also discussed. In 
this Section, we will briefly present these results. 


4.1. Ground states and PECs around the minima 


We studied the ground states of nuclei 74°Cm, ?48Cf, 
250m and 7°2No by using the MDC-RMF model with the 
relativistic density functional DD-PC1 [269]. These four 
nuclei are all well deformed with $29 ~ 0.3 and B49 ~ 0.1. 
But their 822 and (39 values are zero. 


In Fig. 12, the PECs [E = E(32)] are shown for 74°Cm, 
248Cf 250Fm and 2°?No. In the PECs, there exist clear 
potential pockets with the depth more than 0.3 MeV for 
248 Cf and ?50Fm. For 746Cm and 752No, only a shallow po- 
tential pocket appears along the 832 shape degree of free- 
dom. For 74°Cm, the ground state deformation 832 ~ 0.02. 
The PEC is rather flat around the minimum. We have de- 
fined a quantity Egepth, i-e., the energy difference between 
the ground state and the point with 832 = 0, to measure 
the energy gain due to the 632 distortion. For ?4°Cm, 
Eaepth = 34 keV, indicating a very small Y32 correlation. 
For 748Cf, 2°°Fm and 2°?No, the potential minima ap- 
pear (32 ~ 0.037, 0.034 and 0.025 and the energy gain 
Eaepth = 0.351, 0.328 and 0.104 MeV. One can conclude 
that strong Y32 correlations exist in these nuclei. 


Both {32 and Egeptn reach maximal values at 748Cf in 
these four nuclei along the N = 150 isotonic chain, hinting 
that the Yz correlation is the strongest in 74°Cf [125]. 
This conclusion is consistent with that given in Refs. [256, 
267] and with the experimental fact that the 27 state in 
248Cf is the lowest one among those in these four nuclei 
[266]. 
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4.2. Single-particle level structure 


The non-axial octupole Y32 correlation originates from the 
coupling between a pair of single-particle orbitals with 
Aj = Al = 3 and AK = 2. If such a pair of orbitals 
are nearly degenerate and near the Fermi level, one can 
expect a strong Yz effect. 


The single-particle states of ?45Cf around the Fermi lev- 
els are shown in Figs. 13 (for protons) and 14 (for neu- 
trons). In Fig. 13, it can be seen that the two proton 
levels 72 f7/2 and 71%13/2 are close to each other and this 
quasi degeneracy may result in octupole effects. When (29 
increasing from 0 to 0.3, the two proton orbitals [521]3/2 
and [633]7/2 keep nearly degenerate. Since the quantum 
numbers of these two orbitals fulfill the conditions for Y32 
correlations, i.e., Aj = Al = 3 and AK = 2, the cou- 
pling between them leads to non-axial octupole Y32 cor- 
relations. With non-axial octupole deformation (32 in- 
creasing from zero, an energy gap gradually develops at 
Z = 98. In Fig. 14, one finds that the two neutron orbitals 
[622]5/2 (originating from 2g9/2) and [734]9/2 (originating 
from 1j15/2) are also nearly degenerate and just lie around 
the Fermi surface. This results in an energy gap around 
N = 150 with (3 increasing. Thus it was concluded that 
the Y32 correlations are both from protons and from neu- 
trons in these even-even N = 150 isotones and such cor- 
relations are the most pronounced in 748Cf [125]. 


From the above discussions and those made in Sec- 
tion 3.3.3, it can be seen that a proper description of the 
single-particle structure is crucial for the appearance of 
the third barrier, the 632 minimum, et al. Therefore a re- 
lated issue becomes very crucial: How well can covariant 
density functionals (or non-relativistic density function- 
als) describe single-particle level schemes in atomic nuclei? 
Nowadays this has become a hot topic; more details can 
be found in Refs. [270,271] and references therein. 


4.3. Further discussions 


In Ref. [125], the dependence of the above conclusions 
on the functional form and on the effective interaction 
has been studied. Namely, we investigated 74°Cm, ?4°Cf, 
250Fm and 7°2No with effective interactions PC-PK1 [222, 
223], DD-ME1 [166] and DD-ME2 [248]. The results with 
different parameter sets are similar and several general 
conclusions were drawn: (1) Density-dependent function- 
als give stronger Y32 correlations than the nonlinear cou- 
pling ones; (2) The non-axial octupole Y32 effects in ?4°Cm 
are very weak; (3) The Y32 correlations are the strongest 
in 748Cf among these four N = 150 isotones except for 
that DD-ME1 gives a smaller energy gain for 748Cf than 
for 2°°Fm. 


It is worthwhile to mention that the potential pocket in 
the E ~ 632 curve is not quite deep for these four N = 
150 isotones. Therefore the existence of a minimum with 
non-zero $39 in the E ~ (39 curve does not necessarily 
mean that the corresponding nucleus has a static non- 
axial octupole shape. But such a pocket could result in 
Y32 correlations and in lowering the energies of the {32 
vibrations. 
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5. Summary and perspectives 


Many shape degrees of freedom play crucial roles in 
the ground and low-lying states and fission properties of 
atomic nuclei. It is important to include as many shape 
degrees of freedom as possible in a self-consistent model 
for the study of nuclear shapes and potential energy sur- 
faces. This is particularly true for investigating atomic 
nuclei in unknown regions, e.g., exotic nuclei and super- 
heavy nuclei. 

We have developed multidimensionally-constrained co- 
variant density functional theories (MDC-CDFT) [68, 123, 
130,157]. In the MDC-CDFTs all shapes characterized 
by 6), with even u are considered. The covariant den- 
sity functional can be one of the following four forms: 
the meson exchange interaction or point-coupling nucleon 
interaction combined with the non-linear couplings or 
density-dependent couplings. The pairing correlations are 
treated with the BCS approach (MDC-RMF) or the Bo- 
goliubov transformation (MDC-RHB). The MDC-CDFTs 
have been applied to the study of potential energy surfaces 
and fission barriers of actinides [123, 124, 130-132, 158, 
159], the non-axial octupole Y32 correlations in N = 150 
isotones [125], the third minima and triple-humped bar- 
riers in light actinides [133] and shapes of hypernuclei 
[160,161]. In this Review, we presented the formalism 
and some applications on normal nuclei. 

In the study of fission barriers and potential energy sur- 
faces of actinides, we found that besides the reflection 
asymmetric deformation, the non-axial deformation is also 
important for the second fission barriers in actinide nuclei. 
Both the outer and the inner barriers become lower when 
the triaxial distortions are included in the MDC-RMF cal- 
culations. With the inclusion of the triaxial deformation, 
a good agreement between the calculated heights of outer 
fission barrier and the empirical ones is achieved. 

We explored possible existence of the third minima and 
the third barriers in some light actinides and found that 
the appearance of the third minima is dependent on the 
functionals used in the MDC-RMF calculations. The func- 
tional DD-ME2 predicts more prominent third minima 
and barriers than PC-PK1. A detailed analysis of the 
single-nucleon levels around the Fermi surface reveals that 
the formation of the third minimum mainly originates 
from the Z = 90 proton energy gap at 699 ~ 1.5 and 
339 ~ 0.7. 

The non-axial reflection-asymmetric (32 shape in 
246Om, 248Cf, 250Fm and 7°?No were investigated. The 
origin of the non-axial octupole Y32 correlations in these 
nuclei was analyzed by examining the proton and neu- 
tron single-particle orbitals around the Fermi level. The 
Y32 correlations originate from both protons and neutrons 
and such correlations are the most pronounced in 74°Cf 
among these four nuclei. 

At present and in the future, we are carrying out or will 
perform the following investigations: 


e One of our motivations to develop the MDC-CDFTs 
is to study shapes and PES’s of superheavy nuclei. 
We are performing MDC-CDFT calculations for su- 
perheavy nuclei with the axial and reflection symme- 
tries simultaneously broken. 
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e By constraining the shape to an extremely large de- 
formation along a certain shape degrees of freedom, 
one may force some a-nuclei to be in a cluster con- 
figuration, such as the linear chain shape obtained by 
constraining 32 > 0 [272,273], the tetrahedral shape 
by constraining (632 [274-276], the bubble shape by 
constraining the radius, the toroidal shape by con- 
straining 2 < 0, etc. The MDC-CDFTs provide a 
useful tool for the study of such constraint cluster 
structure in atomic nuclei. 


e Based on the MDC-CDFTs with one or more hy- 
peron(s) included, we will study various shape ef- 
fects in light and medium-heavy hypernuclei: (1) Ex- 
otic shapes, e.g., superdeformed and hyperdeformed 
states, the tetrahedron-like state, et al.; (2) The in- 
terplay between localization (clustering) effects and 
spontaneous symmetry breaking of nuclear shapes; 
(3) The effects of hyperon on fission properties; (4) 
The dependence of the shape polarization effect of 
hyperon and the persistence of superdeformed and 
hyperdeformed states on the effective interactions. 


e Most of our investigations presented here have been 
performed for a limited number of nuclei and with 
only several selected effective interactions. Nowa- 
days systematic studies of the global performance 
of many relativistic functionals on the ground state 
properties and beyond-mean-field correlations of nu- 
clei across the nuclear chart have been performed, see, 
e.g., Refs. [277,278]. Such global analysis of the per- 
formance of various functionals on nuclear shapes, po- 
tential energy surfaces and fission barriers will be also 
carried out by using the MDC-CDFTs. 
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